Ultrafast Dynamics of Carrier Multiplication in Quantum Dots 
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A quantum-kinetic approach to the ultrafast dynamics of carrier multiplication in semiconductor 
quantum dots is presented. We investigate the underlying dynamics in the electronic subband 
occupations and the time-resolved optical emission spectrum, focusing on the interplay between 
the light-matter and the Coulomb interaction. We find a transition between qualitatively differing 
behaviors of carrier multiplication, which is controlled by the ratio of the interaction induced time 
scale and the pulse duration of the exciting light pulse. On short time scales, i.e., before intra- 
band relaxation, this opens the possibility of detecting carrier multiplication without refering to 
measurements of (multi-)exciton lifetimes. 

PACS numbers: 78.47.jd, 78.67.Hc, 73.22.Dj 



I. INTRODUCTION 

Typically, in a semiconductor excited by a single pho- 
ton, the conversion efficiency of light energy into electri- 
cal energy is strongly limited by the fact that per photon 
only a single electron-hole pair (or exciton) can be cre- 
ated. For single-junction solar cells, this determines the 
well-known Shockley-Queisser limiti Using multiple ma- 
terials with different band gaps (multi-junction cells), a 
larger part of the solar spectrum can be covered, and 
thus the conversion efficiency can be increased beyond 
this limit. Another way is to use materials which ex- 
hibit the creation of multiple excitons per incident pho- 
ton (multiple-cxciton generation - MEG), thus increas- 
ing the quantum yield for incident photons with energies 
higher than twice the band gap. For this reason, carrier 
multiplication is currently the subject of intensive inves- 
tigation in nanostructures such as graphene structures^ 
and semiconductor quantum dots (QDs)4^— 

Most of the experimental investigations focus on the 
decay dynamics of the carrier population in the con- 
duction band as a measure of MEG: While an exciton 
decays radiatively on a time scale of tens to hundreds 
of picoseconds, multiexcitations can decay much faster 
(on a picosecond time scale) due to additional Coulomb- 
mediated impact ionization and Auger scattering chan- 
nels. Thus, the detection of different decay time con- 
stants in dependence on the incident light frequency is 
used as an indication and quantitative measure of MEG. 
One disadvantage of these approaches is the ambiguity 
accompanying the discrimination of MEG from other ef- 
fects in quantum dots such as charged excitons^— 

Corresponding to experimental investigations, theory 
has focused on calculations of (multi-)exciton decay rates 
which were found to be crucial for MEG and its efficiency 
in concrete material systems 2& The theoretical descrip- 
tions are typically restricted to Bloch equations includ- 
ing phcnomcnological scattering rates as well as pseudo- 
potential approaches without quantum-kinetic time res- 



olution of the carrier multiplication process. The time 
dynamics of absorption bleaching was investigated by 
Shabaev et al. in Ref. H. 

Here, we present a quantum-kinetic approach to car- 
rier multiplication in a semiconductor QD including the 
time dynamics of the relevant Coulomb processes (impact 
ionization and Auger recombination), studied in the elec- 
tronic occupations and the time-resolved optical emission 
spectrum. We propose that distinct quantum-kinetic sig- 
natures of the carrier dynamics due to different dynam- 
ical paths can be used as an experimental verification of 
carrier multiplication via the time-resolved optical emis- 
sion in quantum dot structures: This proposal is, in prin- 
ciple, advantageous compared to measurements of exci- 
ton or multiexciton decay times due to a clear interpre- 
tation of the signal. 



II. MODEL SYSTEM 

The QD model underlying this work consists of eight 
electronic states exhibiting spin degeneracy, i.e., four in 
the conduction (c) and four in the valence band (v), de- 
picted schematically in Fig. [IJa) . The conduction band 
states are labeled by \c\s), while the valence band states 
are labeled by |t>As); here, A denotes the subband index 
(two subband states for each band) and s the spin index. 

The full system Hamiltonian reads 

H = H o + H c i + -ffcoui, (1) 
where Hq represents the free kinetics of the electrons: 
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) and a cXs (a vXs ) are the creation and annihila- 
tion operators of an electron in the conduction (valence) 
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FIG. 1: (a) Quantum dot model in its ground state con- 
sisting of eight electronic states [four in the conduction (c) 
and four in the valence band (v)] exhibiting spin degener- 
acy, (b) Creation of an electron density in the state |c2t) 
by a classical light pulse, (c) Impact ionization (black) and 
Auger recombination (grey) due to energy-conserving con- 
duction band electron-electron scattering, (d) Impact ioniza- 
tion (black) and Auger recombination (grey) due to energy- 
conserving electron-electron scattering in the valence band, 
(e) Simultaneous excitation of two electrons from the valence 
to the conduction band. 



Hamiltonian reads 



,t 



a c2t a 'i>H a eH a clt' 



with v 11 = v? t £nVr 



(4) 

(V AH )*. The Coulomb coupling 

elements V u and V AR are taken to be real and are thus 
equal, V u = V AR = V. 

The theoretical description developed in this work is 
inherently non-Markovian and hence allows to study 
the dynamics including non-energy-conscrving scattering 
processes. This will be utilized later to study the effect of 
off-resonant Coulomb scattering on the efficiency of the 
impact ionization process. 

For our calculations, we assume that the coupling ele- 
ments are determined by an ab initio theory such as in 
Ref. [H, thus focusing on the ultrafast time scale of the 
quantum-kinetic dynamics. The parameters are chosen 
for the PbSe QD of Ref. yielding approximate values 
for the energetic and coupling parameters. 



band and subband A with spin s and energy e^ s , respec- 
tively 

Initially, the systems is excited by a classical light pulse 
close to the transition |«2f) — Hc2f) [Fig. QJb)], creating 
the single-exciton state |c2f v2 J [): This interaction be- 
tween the electrons and the incident light pulse is treated 
semiclassically and in dipole and rotating wave approxi- 
mation: 

fld = M^E(t) al 2f a v2t + h.a. , (3) 

where E(t) is the classical light field at the QD 
location, taken to be a Gaussian pulse E{t) = 
E(t) cos(u L t) with E(t) = E exp{-i 2 /r 2 }, with the 
pulse width r. The pulse area is defined via the enve- 
lope of the pulse, 9 = dt'M^ 2t E(t')/H. Mgfa is the 
dipole matrix element of the considered transition. 

To focus on the underlying dynamics of MEG, the full 
carrier-carrier interaction Hamiltonian is restricted for 
large parts of the article to the impact ionization (II) and 
Auger recombination (AR) processes of conduction band 
electrons, depicted in Fig.njc). The additional consider- 
ation of II and AR of valence band electrons [Fig. [T{d)], 
which occurs on an equal footing, complicates the presen- 
tation of the underlying physics while it does not quali- 
tatively change the observed phenomena. For complete- 
ness, we include calculations of both scattering processes 
at the end of the numerical discussion including a com- 
parison between the system dynamics excluding and in- 
cluding II and AR processes of valence band electrons. 

The impact ionization process annihilates electrons in 
the states |c2f) and and creates electrons in the 

multi-carrier state |clf clj,) [black arrows in Fig. [He)], 
thus describing the transition between the single-exciton 
state |c2f u2f) and the two-exciton state [elf cl\. ulj, 
v2'\). The Auger recombination as the inverse process 
is described by the Hermitian adjoint [grey arrows in 
Fig. [He)]. The corresponding carrier-carrier interaction 



III. DYNAMICAL EQUATIONS 

The quantum-kinetic dynamics is derived within an 
equation of motion approach^ The finite character of 
the four-electron system is used to close the hierarchy of 
expectation values arising from the many-particle inter- 
action Hamiltonian Hq ou \, making the calculated results 
essentially non-perturbative with respect to the Coulomb 
interaction. The restriction to four electrons is a good as- 
sumption if the interaction between the electrons inside 
the quantum dot and the electrons of the surrounding 
substrate is negligibly weak^ Within this assumption, 
all correlations containing ten or more electron operators 
vanish because they describe correlations between more 
than four electrons. 

Since we focus on ultrashort time scales (subpicosec- 
ond), the influence of radiative decay on the carrier 
dynamics, which acts on comparibly long time scales 
(Ti > 1 ns) resulting from the small energy gaps, will 
be neglected throughout this work. However, pure de- 
phasing processes are considered via a phenomcnological 
constant 7pd- 

It turns out that for the system introduced in Sec. HT1 
including the impact ionization and Auger recombination 
of only the conduction band electrons, a reduction to a 
two-electron description yields results close to the full cal- 
culations involving all electron correlations. Therefore, 
to achieve greater clarity, we will use the two-electron 
level for a discussion of the system properties. The fi- 
nal dynamical results presented in the figures, however, 
are calculated using the full set of equations, cf. the ap- 
pendix for the equations. We would like to note that an 
analogous set of dynamical equations was used in Ref. H 
to discuss the time dynamics of absorption bleaching in 
QDs. 

To present a clear depiction of the processes and their 
mutual dependencies, the operator abbreviations used 
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are listed in table fl] 



Abbreviation 



Operator 



/V 
\s 

Pde 



A 



vXs vXs 

t „ 



a clt a cl4- a -ul4- a -u2t 
a elt a elJ. a DlJ. a c2t 



TABLE I: Operator abbreviations: /? is the occupation op- 
erator in subband A with spin s in the conduction (valence) 
band. P2t describes the transition of an electron from state 
\v2f) to |c2f> [cf. Fig.QJb)]. The operator of the double exci- 
tation pde [cf. Fig.QJe)] annihilates two electrons in the states 
|w2'[~) and \vl\) and creates two electrons in the states |clt) 
and |cl^}. The operator of impact ionization pu [cf. Fig.[TJc)] 
describes the correlated transition between the states 
|cl|) and the states |c2|), |clt). 



A. Energetic Structure and absorption of the QD 

Before we investigate the full dynamics, we first char- 
acterize the electronic system via the linear optical ab- 
sorption: 



detuning Ae = Huj — (e^ — e^) = between the incident 
light field and the single-exciton energy (e^ — e^.). For 
non-vanishing and sufficiently large Coulomb coupling, 
the detuning A = (e| t - e^) - [(e^ - e^) + (e^ - e?J] 
between the energy of the single-exciton generated by 
the light field [Fig. [ljb)] (ej^, — e^,) and the energy of the 
two-exciton generated by the double excitation (pde) or 
impact ionization (pu) [Fig. [T|c),(e)] [(e^ — e^) + (e^ — 
e\, )] dominates the root of the denominator and therefore 
the maxima of the absorption. 

The linear absorption spectrum of the system for dif- 
ferent coupling strengths V is shown in Fig. E^a). The 
double peak structure is the consequence of the anticross- 
ing as illustrated in Fig. &b) for different detunings A. 
The states are degenerate for vanishing Coulomb cou- 
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(5) 



Since we are interested in the energetic structure of the 
system under excitation of the transition \v2 t) — > |c2 f), 
we can restrict to the microscopic polarization (p2^) [cf. 
tabic H] , which determines the macroscopic polarization 
P(t) via P(t) ~ M^ 2t (p 2t )(t) + c.c. 

The system of dynamical equations defining the linear 
absorption a(uj) via (p2t) reads: 



dt(p 



2t/ 



7pd 



(P2t) 



^(pde) + ±E(t)M$fa, 



(6) 



dt<PDE> = { r [(e? t _ e 2t) + ( e U - e u)] - ^pd \ We) 



FIG. 2: (a) Linear absorption around the energy of the transi- 
tion \v2\) —¥ |c2f) for A = [energy-conserving impact ioniza- 
tion, cf. Fig. QJc)] and different Coulomb coupling strengths 
V . (b) Linear absorption for V = 20 meV and varying de- 
tuning A and energy detuning Ae, showing an anticrossing 
between the single-exciton |c2f v2 J [) and two-exciton state 
|cltc!4.vi;«2t). 



pling; for non-vanishing coupling, new eigenstates which 
are superpositions of the unperturbed states appear. In 
the following discussion, the system will be resonantly ex- 
cited on the higher energetic branch [cf. Fig.JSJb)] except 
where explicitly noted otherwise. This excitation energy 
is found to be Huj ps (e^ — e^t + V) for V > 3 meV 
in the case of a vanishing detuning A, which describes 
energy-conserving impact ionization as considered here, 

cf. Fig. me). 



-V(p2t)- 



(7) 



The analytical solution of the absorption spectrum can 
be obtained via Fourier transformation, yielding 



B. Structure of the dynamics 



a(ui) ~ 3 



Ae lHj PD /^-[(e^-e^ + ^-e^JJ-i/^ypD 

.(8) 

It can be seen that for a vanishing Coulomb coupling, 
the real part of the denominator vanishes for an energy 



Since the process of impact ionization is the fo- 
cus of this study, we start the discussion of the dy- 
namical structure with the expectation value (pu) = 
(fflcit°ci4. a t;:U. a c2t) [cf. table [I], which is the probability 
amplitude of the impact ionization process (generation of 
two electrons in the conduction band via one electron in 
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the conduction and one in the valence band): 

d t {pu) =(jAw (ii » -tpd)(pii) ~ iM(t)(pDB>-l 
+ ^((/2t/u> - <AVu» 



(9) 



with AwW = Tr 1 ^ + - - e£ f ). Here, the 
product of the dipole coupling element and the slowly 
varying envelope of the electric field is denoted with 
M(t) = E(t)M^ 2 ^. In the following, we write the classi- 
cal light field E(t) as a sum of its positive (e luJLt ) and neg- 
ative (e -1 " 1 '*) frequency parts and abbreviate the prod- 
uct of an expectation value (A) and a rotation of n times 
the laser frequency e ±mULt with 



(A) 



±n 



(A)e 



We then apply a rotating wave approximation by ne- 
glecting fast rotating terms for each dynamical equation. 
For example, in Eq. ©, the term ^M(t)(pr>E)+i docs 
not appear because it oscillates approximately with the 
frequency 2wl, whereas the other terms iM{t)(pa^)-\, 
(fzffii), and have no free oscillations. 

The impact ionization (pn) in Eq. Q is driven (i) via 
the optical field M(t) by the double excitation ampli- 
tude (pde)-i and (ii) via the Coulomb coupling V by 
the density-density correlation (f^fii) of two electrons 
in the states |c2f) and \vli) [cf. Fig. Die)]. The cor- 
relation (/If/iY) describes the simultaneous existence of 
electrons in these two levels which is essential for the scat- 
tering process of impact ionization to occur. Because the 
final states of the electrons undergoing impact ionization 
are |clf) and | cl 4.) , a simultaneous occupation of these 
two states (/i-j-Zi^) blocks the impact ionization process. 
The source term due to the double excitation process 
(pde)-i, depicted in Fig. [IJe), an d mediated by the cou- 
pling M(t), depends on the presence of the light field 
and therefore changes (pn) only while the system is in- 
teracting with the light pulse, in contrast to the Coulomb 
coupling V which acts at all times. This fact will become 
important in the analysis of the occuring system dynam- 
ics. 

The dynamical equations for the double excitation 
(pde)-i and the density-density correlations read 



<Mpde)-i =(iAw (DE) -7pd)(pde)-i 

- l -M(t)( P11 ) + l -v{ P2 ^r n ). l 



(10) 



dtif^fu) =?M(*)9«P*r/i4>-i) - faiM* ( n ) 



dt(f C n f C U )=-rVmpu)), 



(12) 



withAw( DE ) = h- 1 (el f + e^-e'l l -e^ f -nujL). The dou- 
ble excitation (pde)-i is driven by the impact ionization 



(pn) via the light field M(t) and by the term (p2t/ii)-i) 
which describes the excitation of an electron from \v2f) 
to |c2f) correlated with an electron density in the state 
\vll). The density-density correlation (f^fi^) is driven 
by this polarization (P2T-/14.) — i an d the impact ioniza- 
tion (pu), while (fi^fi^) is solely determined by impact 
ionization (pn). 

The system of equations is closed with the following 
two dynamical equations: 

dtipnf 11) -1 =(«Ao; (2t) -7PD)(p2f/i4.)-i + zV<Pde)-i 

- iMmryii) (fnfh))> ( 13 ) 



with Aw^fl = R~ 1 (e| t -e v n -hw L ). 



(14) 



A better understanding of the dynamics can be ob- 
tained with the help of a graphical presentation of the 
coupled expectation values, shown in Fig. [3l We start 
from the initial condition that only the density correla- 
tions in the valence band are equal to one, i.e. (/2+/11) = 
1, while all other density correlations and polarizations 
vanish. The dynamics can essentially be divided into 



j4f ^path 1 




\fhno) ' ((P2 T /r i )) 



[ (pde) ) ' 



FIG. 3: Coupling scheme of the system dynamics given by 
Eqs. (P))- (|14|) . Grey arrows depict electron-light coupling and 
black arrows indicate Coulomb coupling. 

two paths, leading from the optically excited polariza- 
tion (p2t/ii)-i t° the impact ionization amplitude {pn), 
which is directly related to the two-exciton state (/f-j-/^): 
one via the density-density correlation (f^fix) (path 1) 
and one via the correlation describing the double exci- 
tation (pde) (path 2). Path 1 can be decomposed into 
two contributions, each coupling two density-density cor- 
relations via a corresponding polarization ((paf/ij,)— ij 
(pn))- The polarization itself is driven by an inversion 
between the corresponding density-density correlations 

((iVu) ~ (/atAV'^Vu) ~ (/rt/u»- If onl y the hght- 
matter-coupling M(t) is considered, the contribution is 
comprised of the polarization (P2t/ij,)-i and the density- 
density correlations (f^fii) and (f^fii)- Including the 
Coulomb coupling V, the contribution consists of the im- 
pact ionization (pn) and the density-density correlations 
(/ff/ij,) and (fiffi^)- Path 2 links the polarizations 
(p2t/ij,)-i and (pn) via the double excitation (pde)-i- 
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If the typical time scales H/V and r (pulse duration) 
of the two interactions V and M(t) obey the relation 
t < H/V . the dynamics can be divided into sequential in- 
teraction steps: For example, for a short light- matter (r) 
and a comparably long Coulomb interaction time (h/V), 
the dynamics is first restricted solely to the left contribu- 
tion of path 1. For longer times, the right contribution of 
path 1 is activated by the non-vanishing density-density 
correlation (f^fi^) > while path 2 remains dormant due to 
vanishing intermittent polarizations. For a light-matter 
interaction strongly exceeding the Coulomb interaction 
time scale, i.e., t 3> h/V, the dynamics is comprised of 
an interference of both paths. 

A case where path 2 dominates the dynamics while 
path 1 remains dormant is not possible, as can be seen 
in Fig. [3] Let us assume that H/V is negligible compared 
to t. Even though (pde) of path 2 is then driven on a 
much shorter time scale than the density-density correla- 
tion (f%ffi±) of path 1 via the polarization (f>2t/i^)j this 
process is compensated by the slower time scale of the 
second step between (pde) and (pu) as compared to the 
case for the density-density correlation (/IV/ii) • Thus, 
no qualitative difference in the behavior is expected when 
taking the limiting case of r >• H/V, which is just due to 
the fact that the Coulomb interaction acts at all times. 

We will come back to the dynamical structure via the 
two paths in the following sections to explain the numer- 
ical results. 



IV. CARRIER MULTIPLICATION 

We now turn to the central focus of the article, the in- 
vestigation of the ultrafast dynamics of carrier multipli- 
cation. As a measure of the process, we define the carrier 
multiplication CM as the ratio of the total carrier occu- 
pation in the conduction band (f c ) and its light-induced 



part (f c ) 



light 1 



i.e., the conduction band carrier density 



which is excited solely by the light field: 



CM 



(fc) 



(15) 



hi 



The quantity CM is equal to one if no additional carriers 
are created by the Coulomb coupling ((/ c ) = (/c)| light )- 
Its maximal value in the model system treated in this 
paper (cf. Sec. [TTJ) is 2. This results from the fact that the 
electronic excess energy is equal to two times the band 
gap energy. We will now derive the explicit expressions 
for the total carrier occupation in the conduction band 
(f c ), its Coulomb- induced part (/c)| Coul and its light- 
induced part (/c) | 

The total carrier occupation in the conduction band 
reads 



(fc 



- ^2( a l\s u c 

X,s 



> = </ 2 c t > + (/ft) + </u> 



(16) 



The occupation (/£, ) is neither driven by the external 
optical field nor by Coulomb scattering and therefore re- 
mains zero. Considering the dynamical equation for (f c ), 



dt(fc) = ^M(t)9f«p2t>-i) + ^V9f«pn)) 



(17) 



dt(fa)\ 



Mfc) 



two source terms (/c) | light and (/c) | Coul can be identified, 
relating to the light-matter coupling M(t) and the impact 
ionization transition probability, respectively. Thus, the 
carrier multiplication CM can be computed. 

In the following, we analyze carrier multiplication 
for both the energy-conserving and the non-encrgy- 
conserving case numerically, using the parameters in ta- 
bic [TT] if not noted otherwise. 



Parameter 



Value 



V V 
e if! £ U 

e it> e i4. 

V 

7pd 



-400 meV 

meV 
400 meV 
800 meV 
20 meV 
1/(500 fs) 



TABLE II: The electronic energy structure and the Coulomb 
coupling elements are taken from Ref. for a PbSe quantum 
dot of diameter d — 5 nm. This system provides a small 
band gap and allows energy conserving impact ionization and 
Auger recombination processes. 



A. General Properties of our Model 

We first discuss the general relation between the total 
carrier occupation and its light- and Coulomb-induced 
parts in our model system. This will be useful in under- 
standing the dynamical behavior of these quantities in 
the next section. 

The inversion driving the impact ionization [cf. Eq. [9] 
and Auger recombination can be expressed in terms of 
the light- and Coulomb-induced parts of the total carrier 
occupation: 



</2 t /n) - (/it/u> = (fc) 



i light 



2(/c 



ICoul 



(18) 



Here, we employed the dynamical equations (|11I12[) and 
the definitions of (f c )l,. and (f c )L , as well as their 

w w | hght w w I Coul 

initial values. For a vanishing inversion, the ratio of 



(fc 



I Coul 



and (f c 



I light 



is obtained as 



(fc 



I Coul 



= l/2(/ c 



I light 



(19) 



This is equivalent to a carrier multiplication CM = 1.5, 
which is reached in a quasi-stationary limit, where all 
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polarizations (in particular, the impact ionization and 
Auger recombination amplitudes) vanish. 

In our model system, this value is obtained on long 
time scales independent of the exciting pulse parame- 
ters, as will be seen in the next sections. This is due to 
the fact that we do not include relaxation effects such 
as phonon-assisted intraband relaxation, since we want 
to focus on the ultrafast time dynamics of carrier multi- 
plication. Furthermore, we will see that, in contrast to 
first expectations, reaching this quasi-stationary state of 
the system does not solely depend on the damping of the 
participating polarizations by pure-dephasing processes, 
but can be reached on an ultrafast timescale by control 
of the pulse duration r. 

A complementary analysis was performed in Ref. [q| for 
the absorption bleaching, focusing on the relation be- 
tween the Coulomb interaction strength and the relax- 
ation rates of the single-exciton and two-exciton states. 
Here, we focus on the relation of the Coulomb coupling 
strength to the externally controllable pulse length, al- 
lowing an external control in the detection of carrier mul- 
tiplication. Since we do not consider relaxation processes, 
we are thus in the regime where oscillations are observed 
in the absorption bleaching in Ref. [^. 



B. Energy-Conserving Case 

We begin our investigation of the dynamics by con- 
sidering energy-conserving impact ionization, i.e., we as- 
sume the energetic structure given in table [TTl where 
e 2t — e it = e U — e H' ^ G differentiate between two dif- 
ferent regimes: (i) the light-matter coupling occurs on a 
shorter time scale than the Coulomb coupling, i.e., the 
pulse duration r is shorter than the typical Coulomb 
interaction time (H/V), t < H/V and (ii) the light- 
matter coupling strongly exceeds the Coulomb coupling, 
r > H/V. 

(i) T i5 H/V: Figure BJa) shows the conduction band 
carrier occupation (/ c ) and the contribution from the 
light- matter interaction (/c)| light f° r an excitation with 



a 20 fs pulse with a pulse area of 0.17T. The light pulse 
induces a density in the conduction band (/ c ) | 1; ht at 

t ps ps. After the pulse, (/c)| light remains constant 
(dashed line) and is exceeded by the total carrier density 
(solid line) (/ c ), which indicates carrier multiplication. 
(f c ) exhibits clear oscillations above the level of (/c)| light - 
As discussed in Sec. IIIIB1 the build-up of the light- 
and Coulomb-induced carrier density occurs sequentially 
via path 1 depicted in Fig. [3] For the considered pa- 
rameter range, on a short time scale, determined by 
the light-matter interaction r, the light field induces 
the density-density correlation (f^fii) which, on a time 
scale determined by the Coulomb coupling H/V, creates 
the density-density correlation via impact ion- 

ization. During this build-up, impact ionization becomes 
increasingly inhibited because the inversion between the 



density-density correlations (/i+/fi) and (f^fii) tends 
to zero, leading to a slowdown of impact ionization until 
a first maximum of (/ c ) is reached, i.e., a state which 
favors Auger recombination over impact ionization. This 
leads to a subsequent decrease of the total carrier occu- 
pation towards a local temporal minimum until impact 
ionization starts to dominate again. The damping of the 
resulting oscillation of the total carrier occupation (/ c ) is 
determined solely by the pure dephasing of the system. 

The basic physical picture of the oscillation is the an- 
ticrossing described in Sec. IIII Bl The spectrally broad 
pulse simultaneously excites both eigenstates, creating an 
inversion between the unperturbed states, i.e., between 

the correlations (/ 2 c t /u) and (/lf/ij.) t cf - E ^ ©]• The 
limit of CM = 1.5 is reached as soon as the pure dephas- 
ing has damped the participating polarizations. 



=> 3 



9 
O I 

• T— I 

§•1 

o 

o o 



</c> 

(/c) light 




1 . 


It = 20fs| 



1.0 -0.5 0.5 
Time (ps) 



1.0 



3 3 



a 

o i 

• i-H 

aj 

Oh 1 

1 

CJ 
CJ 

o 



(fc) 

(/ c ) light 


1 1 








r = 90fs| 

i i 



■1.0 -0.5 0.5 
Time (ps) 



1.0 



FIG. 4: Dynamics of the conduction band carrier occupation 
(f c ) and the contribution due to the light field (/c)| light for a 
pulsed excitation with a pulse area of 0.17T and pulse length 
of (a) t = 20 fs and (b) r = 90 fs. 



(ii) t ^> H/V: Under excitation with a longer light 
pulse (r = 90 fs), the time scale of the light- matter 
coupling strongly exceeds the Coulomb interaction (r 
00 is. H/V ps 35 fs). This prevents carrier multiplica- 
tion, i.e., the creation of the density-density correlation 
(/it/14.))' to occur as a sequential process of initial light 
excitation and following impact ionization alone. Now, 
the paths depicted in Fig. [3] describing a sequential (path 
1) and a simultaneous (path 2) creation of the density- 
density correlation (/{^/^) contribute both and lead to 
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a concurrent growth of the light-induced (/c)| light and 
Coulomb-induced ( f c ) L, , carrier densities in the con- 

w ' I Coul 

duction band. Figure HJb) shows that carrier multipli- 
cation still occurs, even though no oscillations (no sub- 
sequent impact ionization and Auger recombination) are 
found in the total carrier density (/ c ). This is due to a 
balance between the densities driving the impact ioniza- 
tion, and therefore a vanishing inversion, resulting from 
the interplay of paths 1 and 2. This balance causes the 
carrier multiplication limit of CM = 1.5 to be reached 
nearly simultaneously with the exciting light pulse on a 
much shorter timescale than in the case of short pulse 
excitation. 



C. Non-energy-conserving Case 



When the Coulomb-coupled singlc-exciton and two- 



exciton states are off-resonant, i.e., t% 



2f 



■J-t 



-111 



the influence of impact ionization is strongly suppressed, 
because energy conservation is violated by this scatter- 
ing process. To investigate this quantitatively for our 
Coulomb parameters (cf. table HI|) . the time evolution of 
the conduction band densities (f c ) and (f c ) | light is shown 
in Fig. [3a) for an exemplary detuning A = 100 meV 
between the Coulomb-coupled states and a pulse length 
of 20 fs. The detuning A is achieved by varying the en- 
ergy e"i , while the other parameters remain fixed. The 
system is excited resonantly on the single-exciton energy 



c 2t 



-2f 
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FIG. 5: (a) Time evolution of the total carrier occupation in 
the conduction band (f c ) and its light-induced part (/c)| ligllt 
for a detuning of A = 100 meV between the Coulomb-coupled 
states and an excitation with a 20 fs, 0.17T pulse, (b) Impact 
ionization times, defined as the time of the system to reach a 
carrier multiplication factor of CM = 1.316 (= l+0.5(l-e -1 )) 
over the detuning |A| for an excitation resonant on the singlc- 
exciton energy kujL = (elf — e 2t)- 



sidered, which act on a picosecond time scale in the case 
of phonon-assisted relaxation, the carrier multiplication 
will dramatically decrease below CM = 1.5. 

To investigate the influence of off-resonant Coulomb 
scattering processes on the timescale of carrier multipli- 
cation, we plot the time in which a value CM = 1.316, 
corresponding to the exponential growth value of the CM 
maximum, is reached over the detuning A in Fig. [5jb). 
Carrier multiplication times of a few picoseconds, which 
lie in the range of typical phonon-assisted relaxation 
times' occur for a detuning A > lOOmeV. 

Thus, for sufficient detuning between the Coulomb- 
coupled states, impact ionization is expected to be neg- 
ligible. This result can also be used to justify neglect- 
ing impact ionization scattering in the valence band in 
material systems which show a large difference in the 
effective masses in the conduction and valence band. 
Nevertheless, it should be noted that this result is cru- 
cially dependent on the Coulomb coupling strength V. 
A stronger coupling would raise the importance of non- 
energy-conserving impact ionization for a given detun- 
ing A. 



D. Inclusion of valence band impact ionization 

To include the influence of impact ionization and Auger 
recombination in the valence band [Fig. QJd)] in our dy- 
namics, we generalize the Coulomb Hamiltonian Hc ou \: 



-ffcoui =V U a\^a\, , a.,, , a „ + + V AIX a' + al, , a„, ,a 



AR J J 



*clt cl-l- c2t 



c2t"uU"'cl^ u 'clT 



with the coupling elements V u = V%^$ = (V AR )* . 
This introduces additional valence band scattering cor- 
relations as well as correlations linking both bands. 
Figure [5] shows the corresponding dynamics with 
) and without ((/ c ) „) impact ionization of va- 



\c+v 



((fc 

lence band electrons for a short pulse excitation. The 
system is again excited resonant on the higher ener- 
getic branch, and the conduction and valence band im- 
pact ionization coupling elements are taken to be equal, 
V = V = 20 meV. The occuring oscillation describing 
impact ionization and Auger recombination has a higher 
frequency for the inclusion of both processes which corre- 
sponds to a larger effective Coulomb coupling V e s(V, V). 
For the long pulse excitation, no significant differences 
are found. 



The Coulomb-induced conduction band density 
(f c ) I Co j, and thus the process of carrier multiplication, is 
negligibly weak on an ultrafast subpicosecond time scale. 
For very long time scales, the CM again approaches 1.5, 
as discussed in Sec. IIV A[ assuming that no radiative 
recombination and phonon-assisted intraband relaxation 
are considered. However, when these processes are con- 



V. TIME-RESOLVED OPTICAL EMISSION 

In a next step, we investigate the obtained dynamical 
results in the time-resolved optical emission of the quan- 
tum dot during the impact ionization process, opening 
an alternative approach to measure MEG. In order to de- 
scribe this quantum-optical signal, we consider the fully 
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-0.25 0.25 
Time (ps) 



FIG. 6: Time evolution of the total carrier occupation in 
the conduction band with ((fo)\ c+v ) and without ((/ c )| c ) 
impact ionization processes in the valence band. The system 
is again excited on the higher energetic branch. Differences 
between the oscillation frequencies and the mean values of the 
total occupation (/ c ) result from a changed effective Coulomb 
coupling constant. 



quantized light-matter interaction Hamiltonian 



QO 



2t a c2t a u2t C k 



h.a. 



(20) 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 



Time (ps) 
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FIG. 7: (Color online) Time-resolved optical emission spec- 
tra around the energy of the single-exciton e^- = (e^t — e 2t) 
for an excitation with a pulse area of 0.17T and pulse length 
(a) t = 20 fs, (b) t = 40 fs, and (c),(d) r = 90 fs. (a)-(c) 
The time uncertainty of the detector is set to At = 20 fs. A 
transition from oscillating to non-oscillating behavior occurs, 
corresponding to the carrier dynamics discussed earlier, (d) 
For a larger time uncertainty At = 20 fs, the better energetic 
resolution reveals the substructure of the fluorescence corre- 
sponding to the observed anticrossing in the linear absorption. 



with the electron-photon coupling element M. 



2f 



2e V 2f2t 



in dipole approximation; is the pho- 
tonic annihilation (creation) operator in mode k. Analo- 
gous to the consideration of the semiclassical light-matter 
interaction, a rotating-wave approximation is applied on 
the frequency of the external excitation. The hierarchy of 
dynamical equations stemming from the electron-photon 
interaction is truncated on the two-photon level, i.e., two- 
photon assisted electronic correlations are neglected. 

The time-resolved emission spectrum at an observation 
frequency Wk a at a distance z to the detector is given byi2 



ki.k 2 



xe 



-[(^-k^cAt/V^} 2 e -[{k 2 -k s )cAt/V2} 2 



(21) 



with the time uncertainty of the detector At, the quan- 
tization volume Q, the speed of light c, and the vacuum 
dielectric constant eo. The angular frequencies u>i and 
the momenta ki are related via the linear dispersion re- 
lation LUi = cki . The time uncertainty At has to be equal 
to or smaller than the pulse length r in order to resolve 
the time evolution of the exciting light pulse. 

In Figs. Eta-c), the resulting time-resolved emission 
spectrum is shown for increasing pulse lengths r = 20 
fs, 40 fs, and 90 fs. The time uncertainty At is set to 
20 fs in order to resolve the temporal dynamics. In the 
calculations, the full system including impact ionization 
and Auger recombination in both the conduction and va- 
lence band is considered. For a short pulse excitation 



[Fig. Eta)], oscillations in the emission are found, cor- 
responding to the impact ionization-induced carrier dy- 
namics in Fig. @|a) ■ These oscillations become less pro- 
nounced for an intermediate pulse length [Fig. |7£b)] and 
disappear for an excitation with a long pulse [Fig. Etc)]. 
Thus, the quantum-kinetic signatures, resulting from the 
interplay between the light-matter and the Coulomb in- 
teraction, are recovered directly in the time-resolved op- 
tical emission. 

To recover the energetic structure of the QD, the time 
uncertainty of the detector is increased in Fig.Etd) to al- 
low for a sufficient energetic resolution. The observed 
energetic splitting of the fluorescence line corresponds 
to the anticrossing in the linear absorption (cf. Fig. [2]). 
The earlier onset of the signal at the energetically higher 
branch results from the Rayleigh signal of the exciting 
light pulse. 



VI. CONCLUSION 

A quantum-kinetic approach to the dynamics of car- 
rier multiplication in quantum dots was presented, fo- 
cusing on the interplay between the light-matter and the 
Coulomb interaction on ultrafast timescales, where in- 
traband relaxation is of minor importance. The analysis 
of the underlying dynamics revealed two different paths 
leading to multiple carriers in the conduction band, con- 
trollable via the coupling parameters of the two interac- 
tions: the duration of the exciting light pulse and the 
Coulomb interaction time. 

We observe the occurence of oscillations for certain ex- 
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tcrnally controllable excitation parameters. These oscil- 
lations result from the interplay of impact ionization and 
Auger recombination and could in principle be used to 
detect carrier multiplication before relaxation processes 
occur. The corresponding distinct signatures in the time- 
resolved optical emission spectrum open the possibility 
for an alternative approach to measure multi-exciton gen- 
eration in quantum dots. In particular, they might of- 
fer the possibility to distinguish carrier multiplication in 
quantum dots by Auger-type processes from other pro- 
cesses, such as multiphoton excitation, uncontrolled pho- 
tocharging, and influences of the surface chemistry, which 
can lead to ambiguities in the analysis of excitonic decay 
dynamics .— 

It was also shown that the time dynamics does not 
qualitatively change when neglecting impact ionization 
scattering in the valence band, resulting only in small 
modifications in the carrier occupations. Furthermore, 
the investigation of non-energy-conserving impact ion- 
ization and Auger recombination highlights how crucial 
resonance is for the impact ionization process. 



1. Electronic system 



As in Sec. MI 51 we begin with the dynamical equation 
of the impact ionization amplitude: 

dt{pn) =^0rt + e H _ e H _ £ 2t + «fryPD)(pii) 

- |M(t)(p DE )_x + l jV{f^r n ) 

+ ^«/ 2 yu> - </iVu» 

- ^((/If/if/u) + (/2V1V1V) ( A -!) 

Next, we continue with the two-electron correlations: 

1 

rft(PDE)-l = ^( £ lt + £ H _ e U ~ £ 2t _ fi^L + «fi7PD)(PDE)-l 

- l -M{t){p n ) + jyif&ttu 

+ >((/l C T /l>2t)-l - (/l C ;/l>2t)-l) 



:^(/i C T /i>2t)- 



(A.2) 



d t (/ 2 Vu> = ~M(t) 9«A> 2t >-i) - \V 9f((Ru)) 



(A.3) 



Mfitfu) = *V 3((Pn» 



(A.4) 



rft(/l4.P2t)-l =^( £ 2t ~ £ 2t ~ + ^7Pd)(/hP2T>-1 

+ T^((PDE)-1 - (/2 f PDE)-l) 

+ l j;M(t)((fZ t r n ) - (f% t r u )) (A.5) 
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dt(f^fl) = -\M{t) Z{(f im )-i) \v 3((/ 2 >„» 

(A.6) 

The equations necessary to define the total carrier oc- 
cupation and its parts are: 



d t (p2\)-l =T( e 2t - £ 2t _ + ihypr>)(p2t)-l ( A -7) 



'tV c 2t c 2t 

h 



M(t)((f% t )-(f^)) + -V(pu E }- 1 



Appendix: Dynamical equations 

In the following, we present the dynamical equations 
used in the calculations. We restrict the equations to 
impact ionization and Auger recombination processes in 
the conduction band. 



d t (f^)^--M(t)m P2t )^) 



(A- 



dt(f^) = ~M(t) m P 2 t )-i) - \v 9f«Pn» (A.9) 
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dt(f c lt ) = -rV 9f((pn)) 



d t {f c n ) = 



(A.10) 



(A.ll) 



+ ^((/lt/2 t /u> - </xt/u/n» 

+ i (ej t + - - e| t + z 7PD ft) (/ 2 >n> (A.22) 



^t(/lt/l>2t)-l = 7^(/ 2 >De)- 



+ t(4t - e 2t - hbJ L + ^7 PD )(/ 1 c t / 1 c i p2t)-i (A.12) 



(A.23) 



The higher electronic correlations which are not driven 
are not displayed here. 



Finally, the dynamical equations for all higher correla- 
tions arc given by 



+ T ( e 2t - eat - ^ + ^7PD)(/ 1 c i / 1 >2 t )-i (A.13) 



2. Quantum-optical system 

Next, we show the photon-assisted correlations, again 
restricting to the impact ionization and Auger recombi- 
nation processes in the conduction band. 



2t 2t 



j:m){{fttf2tfu) (/2t/it/u)) 
fr^L + i ^ 7 PD)(/ 1 c t / 1 >2t)-i (A.14) 



d t (c 



*\ L k2 L kl/ 



^Af 2 k t (p 2t c kl ) 



fc"^2t (^2t u k2/ 



-^(?*Wk 2 - ^ki)(c k2 c ki) ( A - 24 ) 



^(/2t/u/u> = ?M(t) S*«/ft/i^>-i) (A.15) *<P2t"ki> = ^KAW+i - ^M(t)</| f c kl ) 



dt(f^f c n f c u ) = rM(i) 9((A c t A>2t)-i) (A.16) 



ft 

+ ^(PDEC kl ) - ^M 2 k t *(/| f ) + ^M k *(/ 2 c t /! t ) 



kl/+l 



+ ^( e 2t _ £ 2t - ^kl + i7PDfr)(P2T c kl) (A. 25) 



dtififfiffn) = i M W »«A c tA>t}-i) (A.17) 



rft(/ 2 >DE)-l = -^M(i)(/ 2 >„) + ^(AV^t)-! 

(A.18) 



^(AVhAV = -j:M(t) 9((A C t A>2t)-l) 



+ t^3«/ 2 >ii» (A.19) 



d t(/2t C kl)+l = ^ M W((P2tC kl ) - (P 2 t C kl)+2) 

~ (4 f }+i + + c L )(/ 2 t ' t c kl ) +1 (A.26) 



rf *(A c t c ki)+i = T^(0((P 2 t C kl>+2 - (PatCfci)) 



h 



+ T-(-^kl + Wi)(/2t C kl)+l 



(A.27) 



dtifi^tfh) = ^((AVi>2 t )-i) (A.20) 



^(AtAVu) = -rM(t) 9((A c t A>2 t )-i) (A.2i) 



d *(P2t C kl)+2 = -^(PDE C kl>+2 

+ ^A/(i)(/ 2 c T c kl ) +1 - ±M(t){f$ t c kl ) +1 
i 



(A.28) 



-(e 2 ' t - e 2t - fiwki + 2cj l + i7p D 7i)(p 2t c kl ) + 2 
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dt(PDEC kl ) = --M(t)( Pu c kl ) +1 (A.29) 

c 

1 

dt(puc kl )+i = -^M(t)(pv E c kl ) (A.30) 
+ ^«/ 2 ViV k i> + i - (/n/u c ki)+i) 

dt(PnCki)+i = ^ M (*)(PD E Cki>+2 (A.31) 

+ ^«A c tAVki>+i - (/2 t /i> k i) + i) 

^(A>2 t c kl ) = ^M(t)((/ 2 t ' t /^c kl ) +1 - (/aVuO+i) 

+ ^(pD ECkl )-^M 2 k f(/ 2 c t / 1 v 

+ ^( £ 2t ^ £ 2t ^ ^kl + i7PD^)(/iV2T c ki> (A.32) 

rf *(PDE C kl)+2 = ^(*)(PlI C kl)+l - ^(/uP2t C kl>+2 

(A.33) 



't(/2t/i> k i) + i = ^W((A>2>c kl } +2 - (A>tc kl )) 

+ ^((p n c kl ) +1 - (p\ l c kl ) +1 ) 

+ ^(-?^ kl + ^L)(/ 2 c t /r ;Ckl >+i (A.34) 



d t(AVu c ki) + i = +^^((Pn c kl) + i - (Pn c kl) + i) 
+ ^(-«Wki + c L )(/ 1 c t / 1 c i c kl ) +1 (A.35) 

^(/2 t AVki)+i = -^ 2 k t (A> 2t >+i 

+ ^W((A>2TC kl } - (A> 2t C k l)+2) 

+ £(-fe*d + ^)(/^ t /r |Ckl ) +a (A.36) 

^(A> t 2t c k 1 )+2 = -^(PDEC k l>+2 (A.37) 
+ ^Ai)((/2tA\c kl >+i - (fitfii^+i) 

+ T( e 2t ~ e 2t _ ^kl + 2w L + i7p D ^)(/ l t '^ 1 .C kl ) +2 
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